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摘 要 : 传统 盲 源 分 离 《blind source separation，BSS) 优化 算法 的 应 用 场合 非常 有 限 ， 而 且 分 离 性 能 不 高 ， 为 此 提 
出 了 一 种 新 的 采用 位 置 混沌 重 构 的 入 侵 杂 草 优化 算法 (invasive weed optimization，IWO)， 并 对 其 在 盲 源 分 离 的 应 
用 进行 了 研究 。 新 算法 在 每 轮 更 新 的 初期 驱动 选 出 的 较 优 个 体 向 此 时 种 群 的 最 优 个 体 做 适当 距离 的 移动 ， 这 样 不 仅 
会 增加 种 群 的 多 样 性 ， 避 免 算法 出 现 早熟 ， 而 且 也 能 够 加 快 收敛 速度 。 育 信号 分 离 仿真 实验 证 实 ， 与 标准 IWO、 粒 
子 群 优化 算法 (particle swarm optimization，PSO) 和 自然 梯度 算法 (natural gradient，NG) 相 比 ， 新 算法 的 性 能 优 
势 明 显 ， 收 敛 速 度 较 快 ， 分 离 精度 较 高 。 
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Application of invasive weed optimization based on location chaos 
reconstruction to blind source separation 
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Abstract: Aiming at the limitation of traditional optimization algorithm for blind source separation (BSS) , this paper 
proposes a new invasive weed optimization (IWO) based on location chaos reconstruction, and explores its application to 
BSS. The new algorithm drivers the superior individuals towards the optimal individual of the population at the early stage 
of each update. It can accelerate the convergence speed, increase the population diversity and avoid the premature 
convergence. Simulation experiments for BSS show that the proposed algorithm is more effective than the basic IWO, 
particle swarm optimization (PSO) and natural gradient (NG) , and has faster convergence speed and higher separation 
accuracy. 
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0 引言 点 用 于 BSS, 如 遗传 算法 (genetic algorithm, GA ) 5.51, PSOIG， 
a 4 和 蚁 群 算法 (ant colony algorithm，ACA) tI 等。 智能 仿生 

盲 源 分 离 (blind source separation，BSS) 是 指 在 源 信 优化 算法 是 通过 模拟 自然 界 生 物种 群 的 行为 特点 、 进 化 规 得 

和 混合 系统 特性 未 知 的 情况 下 ， 仅 由 观测 信和 号 (也 被 称 为 泥 和 生存 方式 而 发 展 起 来 的 一 种 演化 计算 技术 ， 是 人 工 智能 一 
合 信 号 ) 恢复 出 源 信号 各 个 独立 成 分 的 技术 ， 是 当前 信号 处 ”个 新 的 研究 分 支 ， 一 般 不 要 求 判 据 的 凸 性 与 连续 性 ， 甚 至 对 
理 领 域 的 研究 热点 ， 其 研究 成 果 已 广泛 应 用 于 图 像 处 理 、 语 解析 表达 式 都 不 要 求 ， 对 计算 数据 的 不 确定 性 也 有 很 强 的 适 


古 


音 识 别 、 生 物 医学 和 机 械 故 障 检测 等 诸多 场合 0 。 应 能 力 。 
独立 分 量 分 析 (independent component analysis, ICA) 是 将 入 侵 杂 草 优 化 算法 (invasive weed optimization，IWO ) 
组 混合 信号 矢量 还 原 成 一 组 统计 独立 的 矢量 ， 使 得 各 矢量 Mehrabian 等 人 中 提出 ， 模 仿 了 自然 界 杂 草 的 种 子 繁殖 、 
成 分 之 间 的 依赖 性 最 小 的 过 程 。 独 立 分 量 分 析 是 讶 信号 分 离 。 扩散 、 生 长 和 竞争 性 消亡 的 殖民 化 过 程 ， 是 一 种 全 新 的 优化 
的 一 种 重要 技术 。 一般 说 来 , 基于 ICA 的 BSS 有 两 个 重要 组 ， 计算 方法 ， 有 具有 很 强 的 鲁 棒 性 和 自 适 应 性 。 与 其 他 和 群 智能 优 
ll 独立 性 判 据 和 和 优化 算法 。 独 立 性 判 据 用 来 规定 或 化 算法 相 比 ，IWO 算法 在 进化 过 程 中 ， 种 群 的 所 有 个 体 都 有 
描述 分 离 信 号 的 统计 独立 性 ， 而 优化 算法 的 作用 则 是 在 判 据 。 机 会 繁殖 出 下 一 代 个 体 ， 适 应 度 越 高 的 个 体 产 生 的 后 代数 量 
函数 空间 通过 某 种 方法 寻找 到 最 优 解 〈 最 佳 分 离 矩阵 ),， 是 决 。 会 越 多 ， 这 种 机 制 使 得 算法 搜索 能 力 更 强 ， 更 符合 自然 进化 
定 BSS 算法 收敛 速度 和 分 离 精度 高 低 的 关键 性 因素 。 因 此 ， 的 法 则 。 同 时 ， 算 法 不 需要 较 大 的 存储 空间 来 记忆 所 有 个 体 
选择 一 种 恰当 的 优化 算法 是 有 效 实现 BSS 技术 的 重要 保证 。 的 进化 轨迹 ， 空 间 复 杂 度 较 小 。IWO 算法 自 诞 生起 就 引起 了 
传统 BSS 优化 算法 基于 自然 梯度 (natural gradient，NG )， 国际 学 术 界 的 广泛 关注 ， 逐 渐 成 为 解决 实际 工程 优化 问题 的 
不 仅 要 求 判 据 函 数 空间 是 凸 的 、 连 续 可 微 的 ， 而 且 无 法 处 理 种 有 力 工 具 00-1。 
多 峰 的 复杂 优化 问题 ， 所 以 梯度 优化 方法 的 应 用 受到 了 很 大 为 了 克服 BSS 传统 优化 算法 的 局 限 性 和 进一步 发 挥 
的 限制 。 IWO 算法 的 优势 , 本 文 提出 了 一 种 采用 位 置 混沌 重 构 的 入 侵 
为 解决 上 述 问 题 ， 智 能 仿生 优化 算法 在 最 近 几 年 逐渐 被 杂 草 优化 算法 (invasive weed optimization based on location 
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chaos reconstruction, LCR-IWO)，, 然后 将 其 应 用 在 讶 分离 中 。 
1 ” 盲 源 分 离 基 本 算法 
1.1 数学 模型 
言 号 的 混合 过 程 可 以 表示 为 
X(t)= As(1)+n(7) (1) 


其 中 :A4 为 mxn 维 混 合 矩 阵 ，s(D)=[s (0),s (7)…,s, (D] 为 源 信 


号 ，1 为 采样 时 间 ，x()=[%(), 坊 ()…, 太 (人 为 观测 信号 ， 是 


源 s 经 过 A 混合 后 的 信号 数据 ，n(2)=[n(D) (D2, (DJ 


为 噪声 ， 一 般 不 予 考虑 。 
盲 源 分 离 的 目标 是 在 没有 任何 先 验 知识 的 情况 下 ， 仅 仅 
根据 观测 信号 x 计 算出 分 离 矩 阵 W， 进 而 得 到 源 信 号 的 估计 


y(D=[y(D)(0 (OO 。y 也 被 称 为 分 离 信号 ， 这 一 过 程 
可 以 描述 为 


由 


y(OD=Wx(D) (2) 
于 源 信 号 和 传输 信道 特性 都 是 未 知 的 ， 所 以 分 离 信号 
y 的 各 分 量 在 幅度 和 排列 次 序 上 具有 随机 性 ， 这 种 现象 被 称 
为 盲 源 分 离 的 模糊 性 或 不 确定 性 [31。 
1.2 信号 的 预 处 理 
在 一 般 情况 下 ， 观 测 信 号 分 量 间 都 具有 相关 性 ， 所 以 在 
分 离 信 号 之 前 需要 对 信号 数据 进行 预 处 理 ， 提 高 算法 执行 的 
效率 。 信 号 预 处 理 最 常见 的 方法 一 是 使 信号 的 均值 为 零 ， 即 
中 心 化 ， 二 是 预 白 化 或 球面 化 09。 
1) 中 心 化 (centering) 
也 称 均值 化 处 理 , 但 并 不 是 所 有 基于 ICA 的 盲 源 分 离 算 


法 都 必须 要 求 这 一 步 又 ,中 心 化 不 会 影响 到 信号 后 续 的 分 离 ， 
而 且 在 信号 恢复 后 ,还 可 以 再 将 均值 向 量 加 到 恢复 的 信号 上 。 
观测 信号 中 心 化 的 表达 式 为 


X=X-E[X] (3) 

2) 预 白化 (prewritening) 

也 被 称 为 白化 (writening) 或 球面 化 〈sphering)， 目 的 
是 让 信号 各 个 分 量 之 间 不 相关 ， 是 盲 分 离 重 要 的 预 处 理 手段 
之 一 。 白 化 过 程 降 低 了 混合 矩阵 的 维 数 ， 减 少 了 被 估计 的 参 
数 个 数 ， 从 而 也 大 大 降低 了 计算 的 复杂 度 。 对 观测 信号 进行 
预 白 化 的 一 个 常用 方法 是 特征 值 分 解 ( eigenvalue 
decomposition，EVD)， 白 化 的 过 程 实质 为 求 白化 矩阵 Vv。 设 


4  … 轧 为 协 方差 矩阵 E[ xx ] 特征 分 解 后 的 特征 值 ，E 是 
此 特征 值 对 应 的 特征 向 量 组 成 的 矩阵 ， 对 和 角 阵 


D=diag[4 4 … 4%]， 则 白化 矩阵 V 的 计算 表达 式 为 
V=EDYE' (5) 
观测 信号 白化 后 的 输出 v 的 表达 式 为 
V=Vx (4) 
3) 独立 性 判 据 的 选取 
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其 中 : 三 表示 ?的 第 7 阶 累积 量 。 


如 果 3 的 概率 分 布 对 称 ， 那 么 =0， 这 时 负 炳 为 
TO)0) G) 
其 中 
(9) = kunt(») 9) 


所 以 ， 峰 度 可 以 看 做 是 负 灶 的 一 种 特殊 情 
2 入 侵 杂 草 优化 算法 
2.1 算法 原理 

IWO 算法 通过 群体 中 优秀 个 体 来 指导 种 群 的 殖民 化 , 运 
用 正 态 分 布 动态 改变 标准 差 的 扩散 方式 将 产生 的 种 子 散布 在 
父 代 个 体 周围 ， 兼 顾 了 全 局 搜索 和 局 部 搜索 。 然 后 ， 采 用 子 
从 与 父 代 生存 竞争 机 制 ， 较 大 限度 地 保留 了 有 用 信息 。 种 群 
通过 这 种 机 制 不 断 进化 最 终 得 到 最 优 解 。 算 法 考虑 了 群体 的 
多 样 性 和 进化 速度 ， 搜 索 能 力 较 强 ， 能 够 较 大 程度 避免 早熟 


NY 
o 


Py 


的 发 生 。 
入 侵 杂 草 算 法 总 体 可 分 为 : 种 群 初始 化 、 繁 殖 、 空 间 扩 
散 分 布 和 生存 竞争 等 四 个 过 程 外 。 


2.2 数学 描述 
2.2.1 种 群 初 始 化 

确定 初始 杂 草 数量 NM, 和 最 大 杂 草 规模 Nw*、 最 大 迭代 次 
数 ierw 、 问 题 维 数 D 、 繁殖 的 种 子 数 上 限 seed 和 下 限 
差 初 始 值 cw 和 最 终 值 cws 、 杂 草 适 应 度 函数 的 
最 小 值 A 和 最 大 值 。 然后 将 产生 的 种 群 个 体 随机 散布 在 
判 据 函 数 空间 。 
2.2.2 繁殖 

杂 草 在 繁殖 过 程 中 所 生产 的 种 子 数量 与 它 的 适应 度 值 成 
正比 ， 这 样 使 得 适应 度 较 小 的 个 体 只 允许 繁殖 少量 的 子 代 ， 
而 适应 度 较 大 的 个 体 则 可 以 繁殖 较 多 子 代 。 规 定 适 应 度 最 大 
的 杂 草 个 体 产 生 的 种 子 数量 为 seed,。， 适 应 度 最 小 的 杂 草 个 
体 产 生 的 种 子 数量 为 seed,s, ， 而 其 他 个 体 生成 的 种 子 数 与 其 
适应 度 值 服从 向 下 取 整 的 线性 关系 。 杂 草 i 生产 的 种 子 数量 
的 表达 式 为 
seed, -aas 元 人 (Seed — Seed,n, ) 十 Seed Jao 
其 中 : foonow 代表 向 下 取 整 函数 ，f 代表 杂 草 i 的 适应 度 值 。 
图 1 表示 杂 草 i 生产 种 子 的 方式 中。 
2.2.3 扩散 

杂 草 产生 的 种 子 以 正 态 分 布 Y(0.D 的 方式 分 布 在 父 代 个 
体 的 周围 ， 迭 代 次 数 和 标准 差 的 关系 如 下 : 


_ (iter,,, -iter) 
(Cin — Cima) + im GD) 


其 中 : iterw 和 iter 分 别 表示 和 连 代 循环 的 最 大 次 数 和 当前 迭代 
数 ，9 为 非 线性 调和 因子 ， oi 为 标准 差 初始 值 ， cnm 为 标 
准 差 最 终 值 。 从 式 (11) 中 可 以 看 出 ， 标 准 差 o 的 逐渐 减少 使 
父 代 的 搜索 区 域 和 种 群 的 多 样 性 被 逐渐 减少 ， 这 样 会 使 适应 


Seed in 、 标 ; 


ite 


峭 度 和 负 焙 是 信号 处 理 领 域 常见 的 统计 独立 性 度量 方 
法 。 通 常情 况 下 ， 用 四 阶 累积 量 定义 随机 变量 的 峭 度 ; 


kurt(7)=BEL7]-3E 7] (6) 
Gat 确 计 算 ， 可 以 通过 下 式 来 近似 替代 05]: 


vs a (7)+ sg (7)- SA (7) (7) 


高 全 0)+ 


度 值 较 好 的 个 体 聚 集 在 一 起 而 一 些 适 应 度 值 较 差 的 个 体 被 淘 
汰 。 
2.2.4 生存 竞争 
经 过 一 定 迭 代 次 数 后 ， 杂 草 种 群 规模 会 达到 设 定 的 上 限 
wx。 此 时 ， 种 群 进行 生存 竞争 ， 其 规则 是 父 代 和 子 代 所 有 
个 体 按 适应 度 值 高 低 进 行 排序 ， 适 应 度 较 高 的 排序 前 ww 株 
杂 草 个 体 将 保留 下 来 继续 进化 ， 而 其 他 适应 度 低 的 杂 草 个 体 
根据 适 者 生存 的 原则 被 淘汰 。 此 后 ， 杂 草 种 群 规模 会 按 这 种 
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方式 维持 不 变 ， 只 保留 生存 色 


李 著 成 ， 等 : 
E 力 强 的 杂 草 进行 后 续 的 搜索 。 


图 1 杂 草 i 生产 种 子 的 方式 
Fig.1 Seed production procedure of weed i 

2.3 位 置 混沌 重 构 策略 

标准 IWO 算法 虽然 具备 很 多 优点 ， 但 仍然 存在 一 些 缺 
陷 : 对 于 高 维 问题 空间 ， 算 法 搜索 的 时 间 较 长 ， 其 次 ， 算 法 
种 群 的 繁殖 和 生存 竞争 机 制 主要 依据 个 体 适 应 度 的 高 低 来 进 
行 ， 适 应 度 小 的 杂 草 只 允许 产生 少量 种 子 或 直接 被 淘汰 ， 可 
能 会 丢失 一 些 附 近 存 在 全 局 最 优 解 的 个 体 ， 某 程度 上 削弱 了 
杂 草 种 子 全 局 扩散 的 可 能 性 ， 使 算法 容易 持 部 。 

为 此 , 本 文 提出 一 种 基于 位 置 混沌 重 构 的 IWO 算法 , 在 
每 轮 算法 迭代 的 初期 预先 对 选 出 的 优秀 杂 草 个 体 进 行 位 置 重 
构 ， 驱 动 这 些 个 体 向 种 群 此 时 的 最 优 个 体 移 动 ， 使 部 分 杂 草 
个 体 跳 出 局 部 极 值 点 , 快速 引导 整个 种 群 搜索 到 全 局 最 优点 。 
文献 [16] 提 出 的 Firefly Algorithm (FA) 思想 为 位 置 重 构 的 实 
现 提供 了 方向 。 
首先 ， 定 义 杂 草 个 体 i 和 j 之 间 的 吸引 力 ( 其 本 质 可 以 
看 成 为 对 方 杂 草 所 在 田野 对 本 杂 草 的 吸引 力 ) 表达 式 为 


E=&xe” (12) 


其中 :5 为 最 大 吸引 系数 ， 因 为 两 株 杂 草 之 间 的 吸引 力 会 随 
着 二 者 距离 的 增加 逐渐 减弱 ， 吸 引 系数 的 引入 用 来 体现 这 种 
特性 ， 广 为 杂 草 个 体 : 和 7 之 间 的 空间 距离 。 

那么 ， 杂 草 i 向 杂 草 j 移动 的 距离 由 下 式 决定 : 


XkK+l)=x (Kk) +Ex (x (KN (kK) try x(rand—0.5) (13) 


其 中 %、 性 为 个 体 i 和 六 对 应 的 空间 位 置 ，w e[0,1] 为 步 长 
天 子 ， rand s[0.1 为 服从 均匀 分 布 的 随机 因子 。 
一 般 将 由 确定 性 方程 得 到 的 随机 性 的 运动 状态 称 为 混 
沌 ， 自 然 界 和 社会 生活 中 混沌 现象 随处 可 见 ， 是 非 线 性 系统 
中 一 种 较为 普遍 的 现象 ， 行 为 复杂 且 类 似 随 机 。 但 看 似 一 片 
混乱 的 混沌 背后 其 实 存在 着 一 定 规律 。 混 沌 现象 具有 随机 型 、 
规律 性 和 遍历 性 三 个 性 质 ,非常 适合 用 来 研究 全 局 优化 问题 ， 
其 是 遍历 性 ， 可 以 帮助 改进 智能 优化 算法 的 不 足 ， 避 免 其 
早熟 [171。 
因此 ， 为 了 进一步 加 强 种 群 的 多 样 性 ， 考 虑 运用 混沌 理 
论 对 式 (12) (13) 的 权重 进行 混沌 处 理 。 混 沌 模型 的 形式 有 
很 多 种 ， 不 同 的 混沌 模型 对 混沌 优化 过 程 的 影响 不 同 。 最 常 
见 的 是 Logistic 混沌 模型 ， 但 Logistic 遍历 的 不 太 均 匀 性 会 
影响 算法 的 性 能 。 立 方 混沌 映射 要 比 Logistic 映射 具有 更 快 
的 收敛 速度 与 更 好 的 遍历 均匀 性 。 所 以 ， 本 文采 用 立方 混沌 
模型 对 和 ww 进行 构造 : 


yi 


SK+l)=4(5 (Kk) -35 (K) (14) 
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y+1)=4(y (0) -3 (hk) (15) 
在 实际 应 用 中 ， 只 要 立方 映射 的 迭代 初 值 不 为 0， 混 沌 
就 会 发 生 。 
3 ”采用 位 置 混沌 重 构 的 入 侵 杂 草 优 化 盲 分 离 算法 
3.1 算法 流程 图 


图 2 为 基于 位 置 混沌 重 构 的 入 侵 杂 草 优 化 (LCR-IWO) 
的 盲 分 离 算 法 流程 图 。 


7 


读 取 观测 信号 ， 并 对 其 进行 
中 心 化 、 预 白化 处 理 


计算 每 株 杂 草 产 生 的 
种 子 个 数 


按照 正 态 分 布 方式 布置 产生 的 
种 子 ， 然后 种 学 成 长 为 新 杂 章 


对 所 有 杂 草 个 体 
按 适 应 值 由 高 到 低 排 序 


否 
达到 最 大 种 群 个 数 ? 
是 


保留 排序 午前 的 NU 不 
杂 草 个 体 ， 消 除 其 他 个 体 


到 2 基于 LCR-IWO 的 BSS 算法 流程 
Fig.2 Flow chart of BSS based on LCR-IWO 
3.2 ”算法 实现 步骤 
根据 图 2 得 出 算法 的 实现 步 又 
a 读 取 观测 信号 x， 对 其 进行 均值 化 、 白 化 处 理 ， 
b) 初 始 化 参数 N,、 Ns 、 iterww 、 seedww 、 seedun 、 Ginia 、 
,' 、 非 线性 调和 因子 6、h, 、 jw ， 随 机 产生 初始 杂 草 种 


c) 计 算 杂 草 个 体 的 适应 度 值 ( 峭 度 )， 并 记录 最 优 个 体 的 
适应 度 值 (x) ; 


d) 对 新 组 建 的 种 群 个 体 进行 位 置 混沌 重 构 ; 
e) 计 算 每 株 杂 草 产 生 的 种 子 个 数 ; 
人 计算 子 代 个 体 正 态 分 布 的 标准 
并 成 长 为 后 代 杂 草 个 体 ; 
外 计算 父 代 和 子 代 杂 草 个 体 的 适应 度 值 ， 并 按 适 应 度 值 
从 大 到 小 排序 ， 如 果 当 前 种 群 最 优 杂 草 的 f 优 于 公告 板 记 录 


的 f(x)， 则 更 新 f(x) ; 


b) 若 新 种 群 的 杂 草 数量 超过 Ns 时， 保留 适应 度 值 前 
Nuu 株 杂 草 组 成 新 的 种 群 ， 此 后 种 群 的 规模 保持 为 Ns; 


兰 ， 生 成 的 种 子 被 散布 


i) 若 友 代 次 数 达 到 了 结束 条 件 ， 算 法 停止 运行 ， 输 出 结 
果 ， 否 则 转 到 d) 继续 执行 。 
3.3 计算 机 仿真 实验 
为 了 验证 新 算法 的 有 效 性 ， 本 节 通 过 讶 信号 分 离 的 仿真 
实验 , 将 提出 的 新 算法 与 标准 IWO、PSO 和 NG 算法 进行 分 
离 性 能 的 比较 。 
3.3.1 实验 环境 及 参数 取 值 
图 3 为 三 路 具有 3000 个 采样 点 的 源 信号 s : * 为 锯齿 波 
信号 、s; 为 一 段 随 机 的 FSK 信号 、s; 正弦 函数 与 余弦 函数 乘 
积 信号 ， 源 信号 对 应 直方 图 见 图 4。 实 验 随机 生成 的 混合 矩 
阵 4 为 


0.0153 0.8645 0.5361 


0.4104 0.4782 0.6164 
0.1435 0.8632 0.3158 


源 信号 s 被 随机 和 矩阵 4 混合 后 产生 的 观测 信号 x (经 过 
和 白化) 及 其 对 应 直方 图 分 别 见 图 5、6。LCR-IWO 算法 各 参 
数 取 值 情 况 见 表 1。 
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图 3 源 信号 


Fig.3 Source signals 
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图 4 源 信号 直方 图 


Fig.4 Histograms of source signals 
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图 5 混合 信号 


Fig.5 Mixed signals 
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图 6 混合 信号 直方 图 
Fig.6 Histograms of mixed signals 
表 1 LCR-IWO 算法 实验 参数 取 值 

Table ] The parameters used in LCR-IWO 

参数 Nu N iter。 


max 


Seed。 Seed Ginita Gaa O 
_ ”可 产生 的 ”可 产生 的 
、 初始 杂 最 大 杂 ”最 大 和 迭 
意义 最 大 种 子 ”最 小 种 子 
草 规模 草 规模 代 次 数 
数量 数量 


取 值 5 20 100 5 2 0.3 0.0005 3 


标准 差 标准 差 最 非 线 性 调 
初始 值 ” 终 值 ”和 因子 


3.3.2 实验 结果 及 分 析 

7 是 新 算法 在 更 新 16 次 以 后 的 分 离 信 号 (对 应 直方 图 
见 图 8),， 可 以 看 出 源 信号 已 经 被 较 好 地 恢复 。 但 将 图 7 与 图 
3 对 比 发 现 ， 分 离 信 号 各 分 量 与 源 信 号 的 对 应 分 量 相 比 不 仅 
顺序 全 部 发 生 了 改变 ， 而 且 相 位 和 幅度 也 大 部 分 有 所 变化 ， 
这 种 情况 是 由 盲 分离 的 不 确定 性 决定 的 。 

9 是 四 种 优化 算法 的 适应 度 值 曲 线 图 ， 看 出 新 算法 比 
标准 IWO、PSO 和 NG 算法 的 分 离 性 能 要 好 ， 收 敛 速度 和 恢 
复 精度 均 为 最 高 ; 收敛 速度 从 大 到 小 的 排序 为 
speedicrwo >speedwo >speedpso>speedxo ,收敛 时 的 适应 度 值 从 大 到 
小 的 排序 为 fitnessicn wo >fitnessiwo >fitnesspso >fitnessNa 。 

表 2 列举 了 四 种 优化 算法 在 某 些 迭代 时 刻 的 适应 度 值 

(各 算法 收敛 时 的 适应 度 值 被 加 粗 显示 ), 可 知 各 算法 的 收敛 

速度 : speedicrjwo=16、speediwo=27、speedpso=37 ，NG 的 收敛 情况 
未 知 ; 收敛 时 的 适应 度 值 fitnessicrwwo=1.0773 、 fitnesswo=0.8778 、 
fitnesspso=0.8270 ，NG 收敛 情况 未 知 。 所 以 从 表 2 能 够 定量 地 
分 析出 LCR-IWO 的 性 能 明显 高 于 其 他 三 种 算法 。 

10 是 新 算法 得 到 的 分 离 信号 和 源 信号 的 散 点 图 , 从 图 
中 同样 可 以 直观 地 看 出 分 离 信 号 各 分 量 ， 与 源 信号 的 对 应 分 
量 相 比较 ， 顺 序 、 相 位 和 幅度 大 部 分 发 生 了 变化 。 
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图 7 分 离 信 号 
Fig.7 Separated signals 
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图 8 分 离 信 号 直方 图 


Fig.8 Histograms of separated Signals 


录用 定稿 李 著 成 ， 等 
1.2 | | | 


40 
迭代 次 数 


60 


图 9 四 种 优化 算法 的 适应 度 | 
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Fig.9 Fitness curves obtained by four optimization algorithms 
表 2 四 种 优化 算法 在 某 些 迷 代 时 刻 的 适应 度 值 


Table 2 Fitness versus iteration number of four optimization 


algorithms 
算法 5 次 10 次 16 次 ”20 次 27 次 30 次 ”37 次 70 次。 100 次 
LCR-IWO 0.7399 1.0058 1.0773 1.0773 1.0773 1.0773 1.0773 1.0773 1.0773 
IWO 0.7467 0.7754 0.8494 0.8516 0.8778 0.8778 0.8778 0.8778 0.8778 
PSO 0.3938 0.6258 0.6807 0.6807 0.6807 0.7058 0.8270 0.8270 0.8270 
NG 0.0611 0.0844 0.1250 0.1633 0.2351 0.2579 0.2871 0.4515 0.5632 


种 群 规模 wo 选 的 


图 10 源 信号 与 分 离 信 号 的 散 点 
Fig. 10 Scatter plots of separated signal 


算法 各 参数 的 取 值 可 以 遵循 一 


= 


页 


不 要 太 大 , 否则 反 和 1 


S versus source signals 


定 的 规则 进行 08 I。 


初始 


j 会 影响 进化 速度 且 不 会 


明显 提高 分 离 精度 。 此 外 ，WN, 的 取 值 应 该 结合 最 大 种 群 规模 


Na 一 起 考虑 ， 要 与 该 值 有 一 定员 
对 于 大 多 数 优化 问题 ， 
值 较 大 有 利于 


不 能 


DL Ns 


曾 加 种 群 的 多 检 
效 提高 算法 的 性 能 ， 而 


和 殖 的 利 


h 子 数 上 限 seed 


的 取 值 要 依据 实际 情况 ， 


E 离 ， 以 提高 种 群 的 多 样 性 


N 应 在 Na 的 1/8~1/2 之 间 取 值 。N,,。、 
， 但 仅仅 依靠 增加 N, 值 并 
曾 加 算法 的 运行 时 间 ， 所 
一 般 在 10~50 之 间 即 可 。 繁 


般 设 为 3~5 的 整数 , 下 限 seed 一 般 


设 为 1 或 2。 对 于 标准 差 , 首先 要 满足 初始 值 cw 大 于 最 终 值 


oma ;实现 从 进化 初期 全 
一 般 情况 下 在 满 足 Oinitia > Ofina 成 立 的 基础 


差 对 目 


局 寻 优 到 进 


化 后 


期 局 部 寻 优 的 转变 。 


， 充 分 考虑 标准 
标 函 数 寻 优 的 影响 程度 。 非 线性 调和 因子 6 对 算法 也 
有 很 大 的 影响 ， 研 究 发 现 该 值 的 最 佳 选择 一 般 是 3。 
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(LCR-IWO )， 该 方法 在 每 轮 从 代 的 初 其 
个 体 进 行 位 置 重 构 ， 使 部 分 杂 草 个 体 跳出 局 部 极 值 点 ， 从 而 
引导 整个 种 群 以 较 快 的 速度 找到 全 局 最 优点 。 所 谓 
就 是 驱动 选 出 
离 。 
的 某 些 权 习 
优化 算法 应 上 


结束 语 


本 文 提出 了 


上 


种 基于 位 


混沌 重 构 的 入 侵 杂 草 算法 
期 首先 对 选 出 的 优秀 


立 置 重 构 ， 
的 较 优 个 体 向 当时 的 最 优 个 体 移动 适当 的 距 
为 了 进一步 加 强 种 群 的 多 样 性 ， 对 位 置 移动 模型 
进行 了 混沌 处 理 。 在 此 基础 上 , 为 了 克服 传统 BSS 
范围 有 限 和 性 能 不 高 的 缺点 ,本 文 将 LCR-IWO 


UU 


CY 


然后 ， 


uy 


算法 应 用 到 BSS, 与 基于 标准 IWO、PSO 和 NG 算法 的 对 比 
实验 结果 表明 ， 新 算法 的 不 仅 收敛 速度 较 快 ， 而 且 分 离 精度 


也 较 高 。 
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